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Abstract 

We discuss a mechanism that induces a time-dependent vacuum energy on cosmological 
scales. It is based on the instability induced renormalization triggered by the low energy 
quantum fluctuations in a Universe with a positive cosmological constant. We employ the 
dynamical systems approach to study the qualitative behavior of Friedmann-Robertson- 
Walker cosmologies where the cosmological constant is dynamically evolving according with 
this nonperturbative scaling at low energies. It will be shown that it is possible to realize 
a "two regimes" dark energy phases, where an unstable early phase of power-law evolution 
of the scale factor is followed by an accelerated expansion era at late times. 
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1 Introduction 



The idea that Quantum Gravity effects can be important at astrophysical and cosmological 
distances has recently attracted much attention. In particular the framework of Exact 
Renormalization Group (ERG) approach for quantum gravity [1] has opened the possibility 
of investigating both the ultraviolet (UV) and the infrared (IR) sector of gravity in a 
systematic manner. 

The essential ingredient of this tool is the Effective Average Action Tk[gfj,u], a Wilsonian 
coarse grained free energy dependent on an infrared momentum scale k which defines an 
effective field theory appropriate for the scale k. By construction, when evaluated at tree 
level, Ffc correctly describes all gravitational phenomena, including all loop effects, if the 
typical momentum involved are of the order of k. When applied to the Einstein-Hilbert 
action the ERG yields renormalization group flow equations |2] which have made possible 
detailed investigations of the scaling behavior of Newtons constant at high energies [314T5]. 
The scenario emerging from these studies, first demonstrated by Weinberg [16] in c? = 2 + e 
dimensions, suggests that the theory could be consistently defined in c? = 4 at a nontrivial 
UV fixed point where the dimensionless Newton constant, g{k) = G{k)k'^, does not vanish 
in the /c — )■ oo limit, i.e. g{k — )■ oo) = g*. As a consequence the dimensionful Newton 
constant G{k) is antiscreened at high energies, very much as one would expect based on the 
intuitive picture that the larger is the cloud of virtual particles, the greater is the effective 
mass seen by a distant observer |17] . 

Recent works have included matter fields [18] and have also considered a growing num- 
ber of purely gravitational operators in the action. In particular, truncations involving 
quadratic terms in the curvature have been considered in [Tn[T9] - [2T| . while higher powers 
of the Ricci scalar have been studied in [22|[23]. In all the investigations the UV critical 
surface has turned out to be finite dimensional {dyy = 3), implying that the theory is 
nonperturbatively renormalizable. At the UV fixed point the theory has a behaviour very 
similar to QCD, being weakly coupled at high energies but the running of the dimensionful 
Newton constant G{k) = g{k)/k'^ in the deep ultraviolet region is a power-law, at variance 
with the logarithmic scaling of QCD. 

A weakly coupled gravity at high energies is expected to generate important conse- 
quences in several astrophysical and cosmological contexts [21], and in fact the RG flow 
of the effective average action, obtained by different truncations of theory space, has been 
the basis of various investigations of "RG improved" black hole spacetimes, [251427] and 
Early Universe models [281433] . 

However, the behavior of the theory is more complicated at low energies, corresponding 
at cosmological scales, because the /3-functions of any local operator of the type ^/gR, 
y/gR^, . . . , y/gR"^ are singular in the IR due to the presence of a pole at A{k)/k'^ = 1/2 



2 



being A the cosmological constant. The presence of this pole is signahng that the Einstein- 
Hilbert truncation is no longer a consistent approximation to the full flow equation, and 
most probably a new set of IR- relevant operators is emerging in the /c — ?■ limit. This 
singular behavior in fact appears in the ERG for nearly all cutoff threshold functions in 
the Einstein- Hilbert truncation, and it is caused by the presence of negative eigenvalues in 
the spectrum of the fluctuations determinant of the gravitational degrees of freedom. As 
discussed in [M] the dynamical origin of these strong IR effect is due to the "instability 
driven renormalization" , a phenomenon well known from many other physical systems 
[35H37] . We shall see that the low energy domain of the theory is regulated by an IR fixed 
point which drives the cosmological constant to zero at very late times A(t — > oo) = 0. 
This new-type of "decaying A cosmologies" is therefore quite different from previous models 
where the time-dependent cosmological constant encodes the effect of the matter creation 
process [38] . 

Astrophysical consequences of the possible presence of an IR fixed point for quantum 
gravity appeared in [321I39VH2] where it was shown that a solution of the "cosmic coin- 
cidence problem" arises naturally without the introduction of a "dark energy field". In 
particular in the fixed point regime the vacuum energy density p\ = A/SttG is automat- 
ically adjusted so as to equal the matter energy density, i.e. = = 1/2, and the 
deceleration parameter approaches q = —1/4. Moreover, an analysis of the high-redshift 
SNe la data leads to the conclusion that this infrared fixed point cosmology is in good 
agreement with the observations [H]. Recent works have also considered the possibility 
that the "basin of attraction" of the IR fixed point can act already at galactic scale, thus 
providing an explanation for the galaxy rotation curve without dark matter [M1H31 - H5] . 
although a detailed analysis based on available experimental data is still missing. Cos- 
mologies discussing the complete evolution from the early Universe to the present time 
have been considered in [16] and in [17]. In particular in [17] it has been shown that the 
"RG improved" Einstein equations admit (power-law or exponential) inflationary solutions 
and that the running of the cosmological constant can account for the entire entropy of 
the present universe in the massless sector (see also [H] for a review). 

The main purpose of this paper is to explore the idea that quantum effects could 
dynamically drive the cosmological constant to zero at late times so that A{t — t- oo) = as 
a result of an explicit dynamical mapping k — k{t) of the renormalization group trajectories 
generated by the unstable infrared modes of the gravitational sector. 

This idea will then be studied using the so-called Dynamical Systems Approach (DSA), 
a technique already used in cosmology by Bogoyavlensky [19] and further developed by 
Collins, and Ellis and Wainwright to analyze non trivial cosmologies (i.e. Bianchi models) 
in the context of pure General Relativity [50]. Some work has also been done in the case of 
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(minimally coupled) scalar fields in cosmology [5T1[52] and, more recently, of scalar tensor 
theories of gravity [53], f{R) theories of gravity [531 - I57] and Hofava-Lifschits gravity [58] . 
Studying cosmologies using the DSA has the advantage of offering a relatively simple 
method to obtain particular exact solutions and to obtain a (qualitative) description of the 
global dynamics of these models. 

In our specific context the DSA allows us to prove that the presence of a singular 
behavior of the RG evolution for the cosmological constant in the infrared put strong 
constraints on the possible RG trajectories. In particular we will show that it is possible 
to realize a scenario where the Universe has a transition from an early unstable phase of 
power-law evolution of the type a oct^ with /3 g]0, 1[ to a de-Sitter phase. 

The structure of the paper is the following. In Section 2 the basic mechanism of the 
dynamical suppression of the cosmological constant by the unstable low energy modes 
is discussed. In Section 3 the dynamical system analysis of the resulting cosmologies is 
presented. Section 4 is devoted to the conclusions. 



2 Instability induced renormalization 

It is interesting to review in detail the main arguments suggesting that the cosmological 
constant A must have a nontrivial scaling at cosmological distances due to QG effects |34j . 
As already mentioned, the key mechanism is the so called "instability induced renormal- 
ization [351437] . In order to illustrate this point let us look to 22^symmetric real scalar 
field in a simple truncation: 

r,[0] = jd'x {\d,<pd^<p + ^ m\k) <p' + ^ X{k) 

In a momentum representation we have 

^?^ = ^=/ + -^(A:) + A(A:)0^ (2) 

so that r^^"* is positive if > 0; but when < it can become negative for 0^ small 
enough. Of course, the negative eigenvalue for = 0, for example, indicates that the 
fluctuations are unstable, and the non-linear evolution of this instability is a "condensation" 
which shifts the field from the "false vacuum" to the true one, the phenomenon that 
produces the instability induced renormalization. In particular, the /3-functions, obtained 
by p-integrals over (powers of) the propagator p"^ + m?{k) + k"^ are regular in the 
symmetric phase {m? > 0) but there is a pole at = —m{k)'^ — provided k"^ is small 
enough in the broken phase {m? < 0). For k'^ \ \m{kY\ the /3-functions become large and 
there the instability induced renormalization occurs. In a reliable truncation, a physically 
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realistic RG trajectory in the spontaneously broken regime will not hit the singularity at 
k"^ = \m{ky\, but rather make m(A;) run in such a way that |m(/c)^| is always smaller than 
k^. This requires that 



In order to avoid the singularity, a mass renormalization is necessary in order to evolve a 
double-well shaped symmetry breaking classical potential into an effective potential which 
is convex and has a ffat bottom, as it emerges from analytical and numerical calculations 
[37] . However, the truncation implied in ([1]) is not enough to describe the broken phase, 
because its RG trajectories terminate at a finite scale fcterm with k"^^^^ = |m(/ctcrm)^| at 
which the /3-functions diverge. Instead, if one allows for an arbitrary running potential 
Uk{4>), containing infinitely many couplings, all trajectories can be continued to /c = 0, and 
for /c — one finds indeed the quadratic mass renormalization ([3]) as discussed in [37] . 

In the case of gravity we can consider a family of "off-shell", spherically symmetric 
backgrounds labeled by the radius of the sphere 0, in order to disentangle the contribu- 
tions from the two invariants / d'^Xy/g' oc (j)'^ and / d'^Xy/g' i? oc 0^ to the Einstein-Hilbert 
flow. It is then convenient to decompose the fluctuation /i^^ on the sphere into irreducible 
components [1] and to expand the irreducible pieces in terms of the corresponding spherical 
harmonics. For h^^ in the transverse-traceless (TT) sector, the operator T^^^ equals, up to 
a positive constant. 



where = g'^'^ V^Vi^ is the covariant Laplacian acting on TT tensors. The spectrum of 
— V^, denoted {p^}, is discrete and positive. Clearly (j4]) is a positive operator if the cosmo- 
logical constant is negative. In this case there are only stable, bounded oscillations, leading 
to a mild fluctuation induced renormalization. The situation is very different for A > 
where, for k"^ sufficiently small, (jl]) has negative eigenvalues, i. e. unstable eigenmodes. 

A consistent calculation including all the components of the metric fluctuation h^^, 
explicitly illustrates this scenario. Following jlO], the /3-functions for the dimensionless 
Newton constant g{k) = k'^~'^G{k) and the dimensionless cosmological constant A(A;) = 
k~'^A{k) can be introduced 



m{kY oc k"^. 



(3) 



V' + 8(/)-2 + fc2-2A(A;), 



(4) 



hg = /3g{g, \) = [d-2 + r]N]g, 



(5) 



dA = /3x{g,X), 



(6) 



where 



/3a 



(2-r/^)A + 4(47r) 



d vd{d+l) 
' [ 4 



(1-2A) 



d g 



r(r2 + l-f) 
r(n + 1) 



(7) 
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being rjisi = —G{k) ^dtG{k) the anomalous dimension. Its explicit expression reads 



24 ' ' 6 



Tin 



Tin + l) 



where n > 1 is an integer related with the regulator and d is the dimension of the spacetime 

Clearly, the allowed part of the (^-A-plane (A < 1/2) in ((Tj) and (jH]), corresponds to 
the situation k"^ > 2A{k) where the singularity is avoided thanks to the large regulator 
mass. When k"^ approaches 2 A (A;) from above the /3-functions become large and strong 
renormalizations set in, driven by the modes that would go unstable at /c^ = 2 A. In this 
respect the situation is completely analogous to the scalar theory discussed above: its 
symmetric phase (m^ > 0) corresponds to gravity with A < 0; in this case all fluctuation 
modes are stable and only small renormalization effects occur. Conversely, in the broken 
phase (m^ < 0) and in gravity with A > 0, there are modes, which are unstable in absence 
of the IR regulator. They lead to strong IR renormalization effects for k^ \ \m{k)'^\ and 
/c^ \ 2A(A;), respectively. 

We are thus led to the conclusion that the instability induced renormalization should 
occur also in this framework as A; — 0, so that to avoid the singularity the cosmological 
constant must run proportional to k'^, 

A{k) = AJ,^ k^ + subleading terms, /c 0, (9) 

with the constant < 1/2 being an infrared fixed point of the A-evolution. If the behavior 
is actually realized, the renormalized cosmological constant observed at very large 
(cosmological) distances, A{k — > 0), vanishes regardless of its bare value. It is important to 
stress that recent investigations based on a conformal reduction of Einstein- Gravity have 
actually found a new IR fixed point which could represent the counterpart, in the reduced 
theory, of the physical IR fixed point present in the full theory [59] . 

As in the case of the scalar field, the presence of an IR pole is signaling that the Einstein- 
Hilbert truncation is no longer a consistent approximation to the full flow equation near 
the IR singularity, and most probably, a new set of IR-relevant operators is emerging ad 
/c — 7- 0. Although we do not have an explicit solution for the RG flow near the unstable 
phase in the case of gravitjil], near the IR fixed point we can always write 

^ = xf + he', k^o, (10) 

being 6 a critical exponent and /12 a constant related with the eigenvalue of the stability 
matrix around the IR fixed point. Its precise value cannot be determined within the 



^Even for the simple scalar theory the non-perturbative investigation of the IR instability requires the 
use of special numerical techniques to correctly resolve the singularity, see [37] for details. 



6 



linearized theory, but we shall see that our analysis will not depend on the actual value of 
/i2- For ^ > the IR fixed point is attractive, while for 6' < is repulsive. In this latter 
case the IR fixed point is a high-temperature fixed point where /3-functions are suppressed 
as — )■ and the fiows stops before reaching A^^. 

In order to close the system we must map the RG fiow onto the Universe dynamics 
so that k — k{t). Clearly, the Hubble scale H{t) would be a natural choice for the 
infrared cutoff k because in cosmology the Hubble length 1/H[t) measures the size of 
the "Einstein elevator" outside which curvature effects become appreciable [IHIIIT], and 
therefore k ~ H{t). For actual calculations we will set 

k{t)=iH{t), (11) 

being ^ a positive number which we expect should be of the order of unity. We are thus 
led to the following A(i/) dependence: 

A = H^^^ X, + H'^Ch2 = AH^ + B H'^, a = 2 + 29. (12) 

Here < a < 2 correspond with a repulsive IR fixed point, while a > 2 represents an 
attractive IR fixed point. The parameters A = ^^A* and i? = /i2 are introduced for 
notation simplicity instead of ^ and /i2 in the following discussion. Remarkably, we shall 
see that our conclusions will be independent on the actual value of B, while the parameters 
A and a will be the only free parameters of our theory. 

3 Dynamical system analysis 
3.1 Basic equations 

From the previous discussion it is then clear that in this framework the cosmological con- 
stant A is promoted to the status of dynamical variable by Eq. f|T2|) . so that A = A{t). 

Let us then specialize f?^,^ to describe a generic Friedman-Robertson- Walker metric with 
scale factor a{t), where we take T^'' = dia.g[p, —p, —p, —p] to be the energy momentum 
tensor of an ideal fiuid with equation of state p = wp being w > a. constant!^. Then the 
quantum gravity "improved" Einstein equations reduce to the following set of cosmological 

^ The dynamical systems analysis that will follow can be easily generalized to the case of non constant 
w by, for example, adding a new variable corresponding to the pressure. However, this information does 
not add any further understanding of the relation between RG flow and dark energy and we will therefore 
limit ourselves to the case of a single fluid with constant barotropic factor w. 
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equations: 



SnGp 



3 ' 



--7rGp{3w 



P 



-3 {1 + w) Hp-^ 

Gr 



A 



(13) 
(14) 

(15) 



SttG" 

where H = a/a and k = —1, 0, 1 with the usual meaning. It is difficult to determine the 
running of G near the IRFP since the /3-functions are singular near /c — 0, as discussed 
before. However one expects that if an IRFP is also present in the running of g so that 
g{k — 0) = the dimensionful coupling constant G{k) grows without bound as A; — 0, 
so that G{t — 7- oo) = oo at late times. 

In the following we will assume that energy momentum tensor of matter field is con- 
served so that we are led to the following "consistency condition" 



G 



A 



(16) 



and the dynamical evolution of G is fixed by this request as a function of A(t). We shall 
see that the late time behavior of G obtained by the assumption ( IT6l) is actually consistent 
with the possibility that running of g in the infrared is determined by an IRFP. 
By using Eq. (fT2|) in the Friedmann equation we find 



3 K 



{3- A) = BH"" + 8nG p, 



[17) 



which resembles the standard GR one. Assuming if 7^ and A 7^ 3 , we define dimension- 
less variables 



X 



8nGp 



y 



K 



3 K 



{3-A)H^' " (3 -A)' " {3-A)H^a^ 

which will be considered to be functions of the logarithmic time M = ln(a/ao) where is 
the value of the scale factor at some reference time. With some algebra the cosmological 
equations can be written as the first order autonomous system 



X 



'iA + 3) -{A-3)w 



X 



2A 

+ \aiA-3) y' + hA-3) (4 
3 o 



a (1 + 3w)) xy, 



A 
3" 



(2 



a 



(1 + 3-0;) X 



K' 



K 



1-^ 

3 



[1 + 3w) X 



2 
2 



A 

1 

3 



{A -3) {l + 3w) x^ - \ [a{A - 3) + 2 A] y 

(19) 
(20) 

(21) 



y 



Ia 

3 
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with the constraint 



1 + K 



X 



0. 



(22) 



In the above equations the prime stands for the derivative with respect to A/". 

The f l22|) allows to eliminate one of the equations of the ( 1T9H2T|) and obtain a two- 
dimensional phase space. We choose here to eliminate ( 1211) to obtain: 



X 



2A 



-{A + 3)-{A-3)w 



X - - {A-3) {l + 3w) x"^ 



A la{A -3) + 2A]y + l-a{A-3)y' + UA-3) {4-a (1 + 3w)) xy, (23) 
3 3 b 



y 





1 

1 + K 



f ) (2-«) 



{1 + 3w) X 



y 



y, 



X 



y- 



(24) 
(25) 



Note that if y = 0, the above system implies y' = and the x axis is an invariant sub mani- 
fold for the phase space. This means that if the initial condition for the cosmological model 
is ?/ 7^ a general orbit can approach y = only asymptotically. As a consequence, there 
is no orbit that crosses the x axis and no global attractor can exist in the phase space. 



3.2 Finite Analysis 

Setting x' = 0, y' = we obtain the three fixed points in Table [H 

Two of these points {A, C) do not depend on the values of the parameters, but one 
(B) has the x coordinate which is a function of A and the barotropic factor w. This fact 
influences also the value of K i.e. the sign of the spatial curvature index k. Merging will 

occur between A and B for A = and between B and C for A = 0. The first value 

1 + w 

for A represents a bifurcation for the dynamical system and we will not treat it here, the 
second correspond to the case in which A does not contain any quadratic term as a function 
of H. 

The solutions associated with the above fixed points can be obtained integrating the 
equation 

ij=(^l-|^ (^-l-^{l + 3w)x + y^H', (26) 

and are listed in Table [H We can see that ^ is a power law solution whose index resembles 
the standard Friedmann solution, but it is modified by the parameter A. It is easy to see 
that this point represents an expansion if A < 3 and a contraction ii A > 3. The point B 
corresponds to a Milne solution and C corresponds to a de Sitter solution. 

The behavior of the energy density and of the gravitational "constant" can be only 
achieved once an assumption on G has been made. Using the ( IT6l) we find that the behavior 
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of G is in general given by combination of powers of t which depend on w and a as it emerges 
from the results depicted inTable [H However, note in particular that for < A < 3 for 
case A and for case B the evolution of G implies a strongly coupled gravity at late times, 
as implied by the IRFP point model [39]. Substitution in the field equations reveals that 
the point A represents a physical solution [i.e. satisfies the (113II15P ] only if = 0, and 
B = while A can take any value. 

The solution for B is physical only for i? = 0, the space curvature is given by 



Co 



A{l + w) ^ 



(27) 



so that this solution is not flat in general, (i.e. n ^ 0). In the special case a = 2, i? is 
not constrained and A ^ A + B in ( 1271) . The solution associated to the point C instead is 
physical only if 

Po = 0, B = i3-A)H',-^. (28) 

It is interesting to note here that, since in the fixed points A and B the parameter B is 
zero, these points represent states for the cosmology that indistinguishable from standard 
general relativity (GR). This, as we will see, will be an interesting feature in the physical 
interpretation of the orbits. The stability of the fixed points can be inferred with the 
standard techniques of the dynamical system analysis and are summarized in Table |2j 

A global picture of the phase space that summarizes the results above can be found in 
Figures [U and O 

There are also other global properties of the phase space that we can deduce from 
equation (126|1 . Specifically the deceleration parameter q can be written in terms of the 
dynamical variables as 

q = -l - ^ = -^[2A + 2{A-3)y + {A-3)x {l + 3w)] . (29) 
This means that the line 

A l + 3w 
1/ = -^ + ^^^, (30) 



A-3 

divides the phase plane in two regions 



%<^ + i^x, (g>0), ^''^ 

^y<:^ + i^x, (g<0), ^''^ 

in which the decelerating factor is only positive or only negative. This allows us to un- 
derstand if a specific orbit includes the transition between a decelerating an accelerating 



10 



Table 1: Coordinates of the fixed points, eigenvalues, and solutions for the system (11911211) . 



Point Coordinates (x, y) 



Eigenvalues 



A 


[1,0] 


B 


(3-A) (1+3 u;) ' ^ 


C 


[0,1] 



|-(A-3) (a-2) (l + tf;),l + 3tf;-A (l+w)| 
{a(l + ?,w) - 2A{w + 1) + A, «(1 + 3'u;) - 2A{w + 1) - A, } 
|-2,i(3-A) (a -2) (1 + ^)| 



A = ^{a{l + 3u7) - 2A{w + l))^ - A{'iw + If {a - 2){{A - 3)u7 + A - 1) 



Point 



Scale Factor 



Energy Density 



Git) 



A 



/3 



2 A 

p po[t to) G - 2(A-3).po{l+-)=^ (A_3)2(t_iQ)2 



4A{l + wy 



(3-A) 



i3 



a = ao{t -to) p = po{t - to) 



~3—3w 



G 



Ait- to) 



3w+l 



4:71 poil + 3w) 



A 



A 



it - to)' 



C 



a = ao exp [Hoit — to)] 
iA - 3)H^ + BH^ = 



p = 



G — Gn 



A = 3H^ 
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phase typical of the dark energy. In addition, substituting the coordinates of the fixed 
points, we can check if they represent a decelerating or accelerating solution. For point 
C the deceleration factor is always negative as expected by the nature of the associated 
solution. Point B instead lies always on the line ( 130|1 and this is again expected by the as- 

sociated solution. Finally, point A represents an accelerated expansion if < A < 3, 

1 + w 

a decelerated expansion for A < and a accelerated contraction ii A > 3. 

1 + w 

3.3 Asymptotic Analysis 

Looking at fl25l) it appears clear that the phase plane is not compact and it is possible 
that the dynamical system fl23H24l) has a nontrivial asymptotic structure. Thus the above 
discussion would be incomplete without checking the existence of fixed points at infinity 
(i.e. when the variables x or y diverge) and calculating their stability. Such points represent 
regimes in which one or more of the terms in the Friedmann equation (fTTI) become dominant 
and should not be confused with the time asymptotics i.e. t — oo. 

The asymptotic analysis can be easily performed by compactifying the phase space using 
the so called Poincare method [BU] • The compactification can be achieved by transforming 
to polar coordinates {r{Af),9{Af)): 

x-^rcos9, y-^r sin 9, (33) 

and substituting r so that the regime r — oo corresponds to TZ ^ 1. Using the 

coordinates ( l33i) and taking the limit 7^ — ?■ 1, the equations ( l23i) can be written as 

71' = — (A — 3)(a;sin2</) + acos20 — a + 4)(3u;cos0 — 2sin0 + COS0), (34) 

9' = -^4^— 4^ [(l + 3w) cos^-2 sin^l (-l + cos2^-sin2^). (35) 
4(p - 1) 

It can be proven that the existence and the stability of the fixed points can be derived 
analyzing the fixed point of the equation for 9 [60]. Setting 9' = 0, we obtain six fixed 



Table 2: Stability of the fixed points for the system ( 1T91I2T1) . When the stability depends 
on the barotropic factor both the characters of the fixed points are shown. See the text for 
more details. 



Point 


A < 

l+w 


a<2 
1+3"' <A<3 

l+w 


A>3 


A < \+^ '" 

l+w 


a>2 
\+_^''" <A<3 

l+w 


A>3 


A 


repeller 


saddle 


attractor 


saddle 


attractor 


saddle 


B 


saddle 


repeller 


repeller 


attractor 


saddle 


saddle 


C 


attractor 


attractor 


saddle 


saddle 


saddle 


attractor 
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(a) Case ^ < 1. 



(b) Case 1 < A < 3. 



— — , 3 ' 

















-4-2 2 4 



(c) Case ^ > 3. 

Figure 1: Phase plots of the system (11 9112 ip for a < 2 and dust (u; = 0). 

points which are given in Table [3l The solution associated with the asymptotic fixed points 
can be derived in the same way of the ones for the finite case and are also shown in Table 
|3l For details on this derivation we refer the reader to the detailed discussion in [51]. The 
solutions associated to the asymptotic fixed points are of two basic types. A first one is 
an exponential growth i.e. a deSitter phase (points A-B) and a second one whose growth 
and decay depends on the values of a and w. In particular they can represent bounces or 
cosmologies in which the deceleration parameter changes sign. 

Finally, the stability of the asymptotic fixed points can be obtained with the standard 
methods of the dynamical system. The results for the first four fixed points (^oo-^^oo) is 
given in Table |H 

The same is not true for the last two fixed points (£^oo, -^oo) whose stability do depend 
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Table 3: Asymptotic fixed points, 9 coordinates and solutions for the system fll9tl2ip . 



Point 



Behavior 



A 
C 

s 

T 





TT 



37r/4 

77r/4 
arctan (i^) 



a 
a 

a 



ao exp[\Cf (t - to)] 
Qq exp 

Co exp 



ao exp 
exp 



(t — to 

(a-l)(f-fo) \ 



2-a 
2-a 



arctan (^-y^) + tt Oq exp 
r = 1 



/ . _ . \r / -3(2-a)w(3t»+2)-3a+2 \ 
''0) y 3(a-2){u;+l)(3u.-l) J 

3(2-a)w(3w+2)-3a+2 y 
3(a-2){«i+l)(3u)-l) y 

4 



(t - to)^ 



9(a-2)«;^+6(a-2)«;-3a+2 



Table 4: Stability of the asymptotic fixed points for the system ( fT9ll2T]) . 



Point 


a<0 


A<3 
< a < 2 


a>2 


a < 


A>3 
0<a<2 


a>2 


A 


attractor 


saddle 


saddle 


repeller 


saddle 


saddle 


Boo 


repeller 


saddle 


saddle 


attractor 


saddle 


saddle 


c 


saddle 


repeller 


saddle 


saddle 


attractor 


saddle 




saddle 


attractor 


saddle 


saddle 


repeller 


saddle 
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(a) Case A <1. 



(b) Case 1< A < 3. 



y 















• 

> ^ 


• 







-4-2 2 4 



(c) Case ^ > 3. 

Figure 2: Phase plots of the system fll9ti2ip for a > 2 and dust {w = 0). 

on w. The character of these points is comphcated by the fact that the value of TZ', whose 
sign is connected to the stability in the radial direction is zero. Fortunately this problem 
can be solved noting that the next to leading term in the full the equation for TZ' is finite 
and different from zero and can give information on the behavior of the function TZ' nearby 
1Z= 1. The stability thus obtained is shown in Tables [5] and |6l 

Now that the asymptotic fixed points and their stability has been determined let us 
look in more detail to their physical interpretation. 

As said before, these points are characterized by one or more variables to become 
infinite. In terms of the definitions f|T8|) . this corresponds to the fact that either the 
quantities in the numerators are infinite or the ones in the denominators approach zero. 
In the first case we are probably seeing some kind of singularity of the model. In the 
second we are seeing a change in sign of the expansion, which in turns corresponds to a 
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Table 5: Stability of the asymptotic fixed point 8^ 







attractor 


repeller 


saddle 






a< 10 A 3 < A < ^"Tn 

Q+IO 






w 


= 


-10<ckOAA>3 


10 < « < A A < 3"-6 

a+lO 


otherwise 






a > A < A < 3 

a+lO 




w 


= 1/3 


a <0AA>3 
a>0A0<A<3 
a <0AA>3 


a< A A < 


otherwise 


w 


= 1 


< a < 5 A 3"-^ < A < 3 

— a— 5 

a > 5 A A < 3 


a<OAA< =^"-6 

a— 5 

a>5AA> =^"-6 

Q — 5 


otherwise 



Table 6: Stability of the asymptotic fixed point J^c 





attractor 


repeller 


saddle 


w = 


a < -10 A A < 3 
10 < a < A 3"-n < ^ < 3 

a+lO 

a > A A > 3 


a< 10 A A > 3""^ 
a > A A < 3""^ 


otherwise 


w = 1/3 


a<0A0<A<3 
a > A A > 3 
a < A 3°-^ < A < 3 

a— 5 


a < A A < 


otherwise 


w = 1 


0<a<5AA>3 
a>5A3<A< 3^ 

a— 5 


< a < 5 A A < 

a— 5 


otherwise 
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maximum or a minimum of the scale factor. It is important to bear in mind that, because 
of our definition of the time variable, when an orbits "reaches" an infinite point the time 
coordinate changes sign and the Universe follows the orbit with a reversed orientation. In 
some sense one can picture this transition as the fact that the asymptotic points are the 
doorway to a mirror phase space in which the orbit orientation and stability of the fixed 
points are reversed. 

In the pure GR framework, where the cosmological equation predict (almost) always a 
monotonic behavior of the scale factor, extrema of the scale factors can occur only at the 
origin of time. We usually talk then of "bounces" or "(re)collapsing universes" depending 
on the sing of the second derivative of the scale factor (or H). In more complex cosmological 
models, like e.g. f{R) gravity [HI], the behavior is more complicated and the scale factor 
can have in principle a series of extrema located at a generic instant. This is the case also 
for the RG cosmologies. However, we will retain the traditional names to indicate such 
features of the scale factor. 

Let us focus, for example, on the conditions to have a bounce i.e. a situation in 
which a = and a > . In terms of the dynamical variables these can be translated in 
the requirement that one has an asymptotic attractor characterized hj H = = a and 
that this attractor lays in the part of the phase space for which — (1 + g) < which is 
determined via ( l30l) . Using this condition it turns out, for example, that for a > only 
I^oo can represent a bounce. The results for the other points can be found in Table [71 

4 Conclusions 

In this paper we have investigated the effect on cosmological scale of a quantum gravity 
related decaying of the cosmological constant. This investigation was performed using the 
Dynamical System Approach that allows a semi-quantitative treatment of the model via 
the construction of a phase space directly connected to the behavior of the equations. The 
phase space is two dimensional and it contains a total of nine fixed point, of which three 
are finite and six are asymptotical. The stability of these points is found to depend on the 
parameters A and a, but not on B which has a marginal position in the entire analysis. 
The determination of the fixed points and their stability allows us to infer the shape of the 
phase space orbits and to deduce the qualitative features of the cosmic histories possible 
in this models. 

Among these, one is particularly relevant in terms of the problem of dark energy domi- 
nation i.e. the presence of a cosmic history characterized by a first phase in which the scale 
factor grows as a power law with exponent included in ]0, 1[ followed by a second phase 
with a faster growth. In terms of the features of the phase space that such scenario can 
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Table 7: Physical interpretation of the Asymptotic Attractors. 



Point 




Bounce 


(Re) Collapsing 


Singularity 




A 






No 


No 


Yes 




Boo 






No 


No 


Yes 




c 






No 


0<a<2, A>3 


Otherwise 










0<a<2, A>3 


No 


Otherwise 


£ 


{w = 


0) 


a 


< -10 A3 < A< 

— a+10 

-10<cK0Ay4>3 


<a < A -I < A < 

S~A 5 

0<A<3, 0<a<2 


Otherwise 




= 0) 


a 


< 10 A 3 < A < 3""^ 
— a+10 

-10<a<0Ayl>3 


< a < ^t^"/ A 1<A<0 

d~A 5 

0<yl<3, 0<a<2 


Otherwise 




W = 


1/3) 




a < A A > 3 


0<a<2A0<A<3 


Otherwise 


J~ 30 


[w = 


1/3) 




a< A A > 3 


0<a<2A0<A<3 


Otherwise 


Soo{w = 


= 1) 


a 


< 10 A 3 < A < 3°-n 
— a+10 

-10<a<0AA>3 


< a < ^+^°^ A 1<A<0 
< A<3,0 <a <2 


Otherwise 


Too{w = 


= 1) 


a 


< -10 A3 < A < ^ 

— a+10 

-10<«<0AA>3 


0<a< 6t^°/ A 1<A<0 

i—A 5 

0<A<3, 0<a<2 


Otherwise 
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be realized if one has an unstable fixed point representing the first phase and an attractor 
representing the seconcifl. Looking at Tables [2] - 16] one can easily see that there is only is 
one set of the parameters {A, a) for which this can be realized: {A < ^ji^, a < 2), for 
which ^ is a repeller. It is reassuring to notice that in this case A is of the order unity 
and, as A = .^^A*, we deduce that also C, is of the order unity being A* < 1/2. 

When these conditions are satisfied there can exist, depending on initial conditions, one 
orbit that starts at A and either go directly to C or bounce at Soo- The first case can be 
seen as a "classical" Friedmann-de Sitter transition. In the second case the cosmological 
evolution could be richer because the solution Soo can indicate, depending on the sign of 
a, a growth whose rate saturates or an expansion phase comparable with a de Sitter one. 
One could interpret this as a "two regimes" dark energy phase. Note that, as we have 
seen, in the fixed points A and B the parameter B is zero these points represent states 
for the cosmology that are indistinguishable from standard general relativity. This means 
that our model in the neighborhood of these points is indistinguishable from the standard 
cosmology. An example of the evolution of some key cosmological quantities compared 
with the standard de Sitter cosmology is given in Figure [3l and in Figure H] we compare 
the additional RG term in f lT7|) with the standard cosmological constant. 

It is important to emphasize that the standard experimental value of Newton's con- 
stant, Gcxp, does not coincide with the value G{k = $,Hq) which is relevant for cosmology 
today. Gexp is measured (today) at fcexp where the length £ = 4oi is a typical 

solar system length scale, say 10^^ m. Thus, in terms of the running Newton constant, 
Gexp = G{k = C,'/isoi), since 4oi ^ -^o- It is only the cosmological quantity G{k = ^H) 
which dynamically evolves in the fixed point regime, not Gexp- This remark entails that the 
dynamical evolution of the cosmological Newton constant in the recent past does not ruin 
the predictions about primordial nucleosynthesis which requires that G{k = ^H^^ci) coin- 
cides with Gcxp rather precisely. In fact, at the time t = tnuci of nucleosynthesis for which 
H = H(t = tnuci), the cosmological Newton constant was indeed G{k = ^-ffnuci) ~ Gexp 
since ctnuci ~ -^rmci ^soi are of the same order of magnitude. 

The phase space structure also sheds light on aspects of the type of RG cosmologies not 
directly related to the problem of dark energy. For example, it is clear that it is possible 
to interpret the transition to a dark era as the beginning of a primordial inflationary 
phase. Our results then imply that there is only one set of values that allows a "graceful 
exit" i.e. the transition from inflation to a Friedmann cosmology. Specifically in the case 
{A < ^ji^, a > 2) we can have this kind of scenario. However, since there is only one 
point which represents accelerating expansion one can see that the RG model can be used 

■^The second point could be unstable as well, but then the set of the initial condition able to realize 
such scenario would be smaller and more difficult to calculate. 
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to represent either the inflationary era or the dark energy era, but not both. 

Finally the analysis of the asymptotic fixed points gives information on the possibility 
of changes in the sign of expansion rates in this type of cosmologies. In general relativity 
such phases are normally associated either to the so-called "bounces" and are associated 
to cyclic Universes (when if > 0) or to recoUapsing universes (when H < 0). In modified 
theories of gravity changes in the sign of H can occur also when the size of the Universe 
is large because the scale factor does not need to be monotonic [61]. In terms of the 
phase space this kind of behavior is associated to the presence of asymptotic attractors 
and the sign of the quantity — (1 + g) as expressed in fl29|l in the fixed point. In particular 
we have that if — (1 + g) > we have a deceleration followed by an acceleration and if 
— (1 + g) < the opposite situation. A quick analysis shows, for example, that only some 
of the asymptotic attractors for very specific values of the parameters can give origin to a 
bounce. 

The previous results point towards a cosmology with interesting features that, in our 
opinion, deserves some more study. In particular one could try to test the transition to 
Dark Age we have found against the Snia data to obtain some more constraints on the free 
parameters. This, and other issues, will be discussed in a following work. 

References 

[1] M. Renter. Nonperturbative evolution equation for quantum gravity. Phys. Rev. D, 
57:971-985, 1998. 

[2] D. Don and R. Percacci. The running gravitational couphngs. Classical and Quantum 
Gravity, 15:3449-3468, November 1998. 

[3] W. Souma. Non- Trivial Ultraviolet Fixed Point in Quantum Gravity. Progress of 
Theoretical Physics, 102:181-195, July 1999. 

[4] O. Lauscher and M. Renter. Ultraviolet fixed point and generalized flow equation of 
quantum gravity. Phys. Rev. D, 65(2):025013 — h, January 2002. 

[5] O. Lauscher and M. Renter. Flow equation of quantum Einstein gravity in a higher- 
derivative truncation. Phys. Rev. D, 66(2):025026-+, July 2002. 

[6] O. Lauscher and M. Renter. Is quantum Einstein gravity nonperturbatively renormal- 
izable? Classical and Quantum Cravity, 19:483-492, February 2002. 

[7] M. Renter and F. Saueressig. Renormalization group flow of quantum gravity in the 
Einstein-Hilbert truncation. Phys. Rev. D, 65(6):065016-+, March 2002. 



20 



M. Reuter and F. Saueressig. A class of nonlocal truncations in quantum Einstein 
gravity and its renormalization group behavior. Phys. Rev. D, 66(12):125001 — h, De- 
cember 2002. 

D. F. Litim. Fixed Points of Quantum Gravity. Physical Review Letters, 
92(20):201301-+, May 2004. 

A. Bonanno and M. Reuter. Proper time flow equation for gravity. Journal of High 
Energy Physics, 2:35 — h, February 2005. 

A. Codello and R. Percacci. Fixed Points of Higher-Derivative Gravity. Physical 
Review Letters, 97(22) :221301-+, December 2006. 

M. Reuter and H. Weyer. The role of background independence for asymptotic safety 
in Quantum Einstein Gravity. General Relativity and Gravitation, 41:983-1011, April 
2009. 

M. Reuter and H. Weyer. Background independence and asymptotic safety in confor- 
mally reduced gravity. Phys. Rev. D, 79(10):105005-+, May 2009. 

P. F. Machado and R. Percacci. Conformally reduced quantum gravity revisited. 
Phys. Rev. D, 80:024020, April 2009. 

E. Manrique and M. Reuter. Bare vs. Effective Fixed Point Action in Asymptotic 
Safety: The Reconstruction Problem. PoS CLAQG08, page 001, 2009. 

S. Weinberg. Ultraviolet divergences in quantum theories of gravitation. In S.W. 
Hawking and W. Israel, editors. General Relativity, an Einstein Gentenary Survey. 
Cambridge University Press, 1979. 

A. M. Polyakov. A Few Projects in String Theory. PUPT-1394, April 1993. 

R. Percacci and D. Perini. Constraints on matter from asymptotic safety. Phys. Rev. 
D, 67(8):081503, April 2003. 

Max Niedermaier and Martin Reuter. The asymptotic safety scenario in quantum 
gravity. Living Reviews in Relativity, 9:5, December 2006. 

D. Benedetti, P. F. Machado, and F. Saueressig. Four-derivative Interactions in 
Asymptotically Safe Gravity. In J. Kowalski-Glikman, R. Durka, & M. Szczachor, 
editor, American Institute of Physics Gonference Series, volume 1196 of American 
Institute of Physics Gonference Series, pages 44-51, December 2009. 



21 



[21] D. Benedetti, P. F. Machado, and F. Saueressig. Taming perturbative divergences in 
asymptotically safe gravity. Nuclear Physics B, 824:168-191, January 2010. 

[22] A. Codello, R. Percacci, and C. Rahmede. Investigating the ultraviolet properties of 
gravity with a Wilsonian renormalization group equation. Annals of Physics, 324:414- 
469, February 2009. 

[23] A. Codello, R. Percacci, and C. Rahmede. Ultraviolet Properties of f(R)-GRAVITY. 

International Journal of Modern Physics A, 23:143-150, 2008. 

[24] A. Bonanno. Astrophysical implications of the Asymptotic Safety Scenario in Quan- 
tum Gravity. ArXiv e-prints, November 2009. 

[25] A. Bonanno and M. Renter. Quantum gravity effects near the null black hole singu- 
larity. Phys. Rev. D, 60(8):084011-+, October 1999. 

[26] A. Bonanno and M. Renter. Renormalization group improved black hole spacetimes. 
Phys. Rev. D, 62(4):043008-+, August 2000. 

[27] A. Bonanno and M. Renter. Spacetime structure of an evaporating black hole in 
quantum gravity. Phys. Rev. D, 73(8):083005-+, April 2006. 

[28] A. Bonanno and M. Renter. Cosmology of the Planck era from a renormalization 
group for quantum gravity. Phys. Rev. D, 65(4):043508 — h, February 2002. 

[29] A. Bonanno, G. Esposito, and C. Rubano. A Class of Renormalization Group In- 
variant Scalar Field Cosmologies. General Relativity and Gravitation, 35:1899-1907, 
November 2003. 

[30] A. Bonanno, G. Esposito, and C. Rubano. Arnowitt Deser Misner gravity with variable 
G and A and fixed-point cosmologies from the renormalization group. Classical and 
Quantum Gravity, 21:5005-5016, November 2004. 

[31] A. Bonanno, G. Esposito, and C. Rubano. Improved Action Functionals in Non- 
Perturbative Quantum Gravity. International Journal of Modern Physics A, 20:2358- 
2363, 2005. 

[32] A. Bonanno, G. Esposito, C. Rubano, and P. Scudellaro. The accelerated expansion 
of the universe as a crossover phenomenon. Classical and Quantum Gravity, 23:3103- 
3110, May 2006. 

[33] A. Bonanno and M. Renter. A cosmology of the Planck era from the renormalization 
group for quantum gravity. In I. Ciufolini, E. Coccia, M. Colpi, V. Gorini, and 
R. Peron, editors. Recent Developments in Gravitational Physics, pages 461 — h, 2006. 



22 



M. Reuter and H. Weyer. Quantum gravity at astrophysical distances? Journal of 
Cosmology and Astro-Particle Physics, 12:1 — h, December 2004. 

J. Alexandre, V. Branchina, and J. Polonyi. Instability induced renormalization. 
Physics Letters B, 445:351-356, January 1999. 

O. Lauscher, M. Reuter, and C. Wetterich. Rotation symmetry breaking condensate 
in a scalar theory. Phys. Rev. D, 62(12) :125021-+, December 2000. 

A. Bonanno and G. Lacagnina. Spontaneous symmetry breaking and proper-time flow 
equations. Nuclear Physics B, 693:36-50, August 2004. 

J. A. S. Lima. Thermodynamics of decaying vacuum cosmologies. Phys. Rev. D, 
54:2571-2577, August 1996. 

A. Bonanno and M. Reuter. Cosmology with self-adjusting vacuum energy density 
from a renormalization group flxed point. Physics Letters B, 527:9-17, February 2002. 

A. Bonanno and M. Reuter. Cosmological Perturbations in Renormalization Group 
Derived Cosmologies. International Journal of Modern Physics D, 13:107-121, 2004. 

E. Bentivegna, A. Bonanno, and M. Reuter. Confronting the IR flxed point cosmology 
with high-redshift observations. Journal of Cosmology and Astro-Particle Physics, 
1:1-+, January 2004. 

M. Reuter and F. Saueressig. Nonlocal quantum gravity and the size of the universe. 
Fortschritte der Physik, 52:650-654, June 2004. 

M. Reuter and H. Weyer. Renormalization group improved gravitational actions: A 
Brans-Dicke approach. Phys. Rev. D, 69(10):104022-+, May 2004. 

M. Reuter and H. Weyer. Running Newton constant, improved gravitational actions, 
and galaxy rotation curves. Phys. Rev. D, 70(12): 124028 — h, December 2004. 

G. Esposito, C. Rubano, and P. Scudellaro. Spherically symmetric ADM gravity with 
variable G and Ac. Classical and Quantum Gravity, 24:6255-6266, December 2007. 

M. Reuter and F. Saueressig. From big bang to asymptotic de Sitter: complete 
cosmologies in a quantum gravity framework. Journal of Cosmology and Astro-Particle 
Physics, 9:12-+, September 2005. 

A. Bonanno and M. Reuter. Entropy signature of the running cosmological constant. 
Journal of Cosmology and Astro -Particle Physics, 8:24 — h, August 2007. 



23 



[48] A. Bonanno and M. Renter. Primordial entropy prodnction and A-driven inflation from 
Quantum Einstein Gravity. Journal of Physics Conference Series, 140(1):012008 — h, 
November 2008. 

[49] O. I. Bogoyavlensky. Methods in the qualitative theory of dynamical systems in astro- 
physics and gas dynamics. Springer Series in Soviet Mathematics, Berlin: Springer, 
1985, 1985. 

[50] G F R Ellis J Wainwright, editor. Dynamical Systems in Cosmology. Cambridge 
University Press, 1997. 

[51] Edmund J. Copeland, Shuntaro Mizuno, and Maryam Sliaeri. Dynamics of a scalar 
field in Robertson- Walker spacetimes. Phys. Rev., D79:103515, 2009. 

[52] Edmund J. Copeland, Andrew R Liddle, and David Wands. Exponential potentials 
and cosmological scaling solutions. Phys. Rev., D57:4686-4690, 1998. 

[53] S. Carloni, S. Capozziello, J. A. Leach, and P. K. S. Dunsby. Cosmological dynamics 
of scalar-tensor gravity. Class. Quant. Crav., 25:035008, 2008. 

[54] Sante Carloni, Peter K. S. Dunsby, Salvatore Capozziello, and Antonio Troisi. Cos- 
mological dynamics of i?" gravity. Class. Quant. Grav., 22:4839-4868, 2005. 

[55] Jannie A. Leach, Sante Carloni, and Peter K. S. Dunsby. Shear dynamics in Bianchi 
I cosmologies with i?"-gravity. Class. Quant. Crav., 23:4915-4937, 2006. 

[56] M. Abdelwahab, S Carloni, and P K. S. Dunsby. Cosmological dynamics of exponential 
gravity. Class. Quant. Crav., 25:135002, 2008. 

[57] S. Carloni, A. Troisi, and P. K. S. Dunsby. Some remarks on the dynamical systems 
approach to fourth order gravity. Cen. Rel. Grav., 41:1757-1776, 2009. 

[58] Sante Carloni, Emilio Elizalde, and Pedro J. Silva. An analysis of the phase space of 
Horava-Lifshitz cosmologies. Class. Quant. Grav., 27:045004, 2010. 

[59] Elisa Manrique, Martin Renter, and Frank Saueressig. Bimetric renormalization group 
flows in quantum einstein gravity, eprint arXiv, 1006:99, Jun 2010. 35 pages, 3 figures. 

[60] S Lefschetz. Differential Equations: Geometric Theory. Dover (New York), 1977. 

[61] Sante Carloni, Peter K.S. Dunsby, and Deon M. Solomons. Bounce conditions in f(R) 
cosmologies. Class. Quant. Grav., 23:1913-1922, 2006. 



24 




Figure 3: Evolution of some key cosmological quantities in RG cosmologies with w = 
A = 1/2, a = 1 (dashed) compared with the ones of GRA (solid). We have chosen as 
initial condition the phase space point {x = l.l,y = 0.1) present in both the phase spaces. 
The index is associated to the value of all the quantities in this point. Upper panel: Semi 
logarithmic plot of H/Hq as function of ln(a/ao). Lower panel: Logarithmic plot of ^nGp 
as a function of ln(a/ao)- The energy density is not plotted because its behavior is exactly 
matched in the two models. 
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Figure 4: Evolution of the term in RG cosmologies with 1/7 = 074 = 1/2,0 = 1, 5 = 2 
(dashed) compared with the A term in GR (solid). We have chosen as initial condition the 
phase space point (x = 1.1, y = 0.1) and the index is associated to the value of all the 
quantities in this point. 
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